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repulsion  u  and  the  bandwidth  r  has  been  previously  reported; 
however,  this  previous  work  has  been  limited  to  metallic  single 
element  conductors.  This  work  extends  the  theory  to  covalently 
bonded  insulators  (and  possibly  semiconductors)  consisting  of 
more  than  one  element. 

Application  of  the  theory  is  made  to  the  0  KW  and  Si  L-.W 
Auger  lineshapes  from  Si02 .  A  high  energy  shoulder  at  511  ev 
in  the  0  KW  lineshape  is  interpreted  as  arising  directly  from 
correlation  effects.  A  peak  at  50  eV  in  the  Si  1*23^  lineshape, 
its  intensity  significantly  underestimated  by  the  previous 
theory,  is  now  accounted  for;  a  peak  at  70  eV  previously  sug¬ 
gested  to  be  a  shake  satellite  is  now  indicated  also  to  arise 
from  correlation  effects.  Both  lineshapes  reveal  a  density  of 
states  primarily  localized  on  a  Si-0  sub-cluster.  The  magnitude 
of  the  hole-hole  repulsion  on  the  sub-cluster  and  between 
neighboring  Si-0  sub-clusters  is  empirically  determined  from  the 
Auger  lineshapes  to  be  ^  11  and  4  eV,  respectively.  The  oxygen 
2p  nonbonding  bandwidth  is  estimated  to  be  v  6  eV,  but  in  light 
of  other  theoretical  and  experimental  results,  our  result  is 
believed  to  be  1-2  eV  too  large.  Reasons  for  our  over-estimate 
are  discussed. 
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FINAL  STATE  CORRELATION  EFFECTS  IN  AUGER  LINESHAPES; 

APPLICATION  TO  Si02* 

David  E.  Ramaker 
Chemistry  Department 
George  Washington  University 
Washington,  D.C.  20052 

Final  state  correlation  effects  in  Auger  lineshapes  are 
considered  within  the  cluster  LCAO-MO-CI  theory  with  a  para¬ 
meterized  Hamiltonian.  A  model  problem  is  solved  analytically 
to  elucidate  the  role  of  final  state  hole-hole  correlation  and 
to  understand  the  localization  of  the  holes  on  rather  small  sub¬ 
clusters  of  the  system.  The  relationship  of  the  correlation 
effects  to  the  relative  magnitudes  of  the  one-center  hole-hole 
repulsion  u  and  the  bandwidth  T  has  been  previously  reported; 
however,  this  previous  work  has  been  limited  to  metallic  single 
element  conductors.  This  work  extends  the  theory  to  covalently 
bonded  insulators  (and  possibly  semiconductors)  consisting  of 
more  than  one  element. 

Application  of  the  theory  is  made  to  the  0  KW  and  Si  L23W 
Auger  lineshapes  from  Si02.  A  high  energy  shoulder  at  511  eV 
in  the  0  KW  lineshape  is  interpreted  as  arising  directly  from 
correlation  effects.  A  peak  at  50  eV  in  the  Si  L23W  lineshape, 

_ 

This  work  was  supported  in  part  by  the  Office  of  Naval  Research. 


its  intensity  significantly  underestimated  by  the  previous 
theory,  is  now  accounted  for;  a  peak  at  70  eV  previously  sug¬ 
gested  to  be  a  shake  satellite  is  now  indicated  also  to  arise 
from  correlation  effects.  Both  lineshapes  reveal  a  density  of 
states  primarily  localized  on  a  Si20  sub-cluster.  The  magnitude 
of  the  hole-hole  repulsion  on  the  sub-cluster  and  between 
neighboring  Si20  sub-clusters  is  empirically  determined  from 
the  Auger  lineshapes  to  be  ~11  and  4  eV,  respectively.  The 
oxygen  2p  nonbonding  bandwidth  is  estimated  to  be  ~6  eV,  but 
in  light  of  other  theoretical  and  experimental  results,  our 
result  is  believed  to  be  1-2  eV  too  large.  Reasons  for  our  over¬ 
estimate  are  discussed. 

I .  INTRODUCTION 

Recently  we  reported  results  of  an  investigation  of  the 
Auger  lineshapes  involving  the  valence  electrons  in  Si02  (1,2). 
in  this  work  we  compared  the  Si  L^L23W,  the  Si  L23W,  and  0  KLL 
experimental  Auger  lineshapes  with  those  lineshapes  predicted 
from  theory.  The  electronic  structure  of  Si02  for  this  work 
was  described  in  terms  of  one  electron  molecular  orbitals  (M.O.) 
on  minimum-sized  clusters  (i.e.,  central  atom  plus  nearest  neigh¬ 
bors  such  as  SiO*  and  Si206+  in  Si02).  The  results  of  this 
comparison  were  gratifying,  with  the  calculated  and  experimental 
Auger  lineshapes  in  generally  good  agreement  in  all  three  line- 
shapes.  Figure  1  summarizes  this  data.  The  only  significant 
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discrepancies  are:  (l)  the  absence  of  the  high  energy  shoulder 
at  511  eV  in  the  0  KW  theoretical  lineshape,  and  (2)  the  sig¬ 
nificant  underestimation  of  the  low  energy  peak  at  50  eV  in  the 
Si  L23W  lineshape. 

That  good  agreement  between  the  cluster  LCAO-MO  theory  and 
experiment  should  be  obtained  is  perhaps  not  surprising.  The 
observed  Auger  lineshapes  in  the  more  covalent  Si02  are  very 
similar  to  the  corresponding  lineshapes  in  more  ionic  solids 
where  the  cluster  LCAO-MO  model  is  expected  to  be  valid.  Thus, 
for  example,  the  S  and  P  L23W  Auger  lineshapes  obtained  from  the 
highly  ionic  Li  salts  are  very  similar  in  all  aspects  to  the  Si 
L23W  lineshape  (3).  Furthermore,  the  0  KW  lineshape  from  the 
more  ionic  MgO  is  essentially  identical  to  the  0  KW  lineshape 
in  Si02  (4).  The  valence  electron  energy  levels  for  Si02  have 
been  probed  by  photoemission,  by  Si  K^  and  L2  3  x-ray  emission, 
and  by  0  Ka  x-ray  emission  (5).  These  data  have  been  well  des- 
cribed  within  the  Si04  cluster  one  electron  LCAO-MO  model  ( l) y 
indeed,  in  our  work,  the  orbital  energies  were  obtained  from  the 
XPS  data  and  the  LCAO-MO  populations  were  adjusted  to  reproduce 
the  XES  data.  Furthermore,  the  atomic  Auger  matrix  elements 
giving  the  Auger  intensity 
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state  holes  in  the  Auger  process)  were  not  evaluated  in  our  theory 
but  were  estimated  from  gas  phase  Auger  data  (2). 

Despite  these  successes  of  the  M.O.  model,  one  can  justifi¬ 
ably  question  the  validity  of  utilizing  the  cluster  one-electron 
M.O.  model  to  study  Auger  transitions  in  covalent  solids.  The 
Auger  process  produces  two  final  state  holes,  whereas  the  x-ray 
emission  and  photoemission  processes  produce  just  one  hole.  In 
general,  these  two  final  state  holes  see  one  another,  a  repulsive 
interaction  U  exist  between  them,  and  their  relative  motion  is 


correlated.  This  interaction  exhibits  itself  in  the  energy  of 
the  Auger  electron. 


E  •  =  E  -E  -E  ' 
cw  c  v  v 


where  the  E's  are  the  one-electron  binding  energies  of  core  and 
valence  electrons,  respectively.  In  Si02»  U  is  substantial,  vary¬ 
ing  from  5  to  20  eV  depending  on  the  central  atom  and  the  localized 
nature  of  the  final  states  involved.  Hole-hole  correlation  may 
also  substantially  alter  the  Auger  intensities.  We  note  that  the 
hole-hole  correlation  effects  on  the  intensities  of  a  single  atom 
are  in  some  sense  already  included  in  our  model  because  the  experi¬ 
mental  atomic  Auger  matrix  elements  were  utilized.  Thus,  only 
the  interatomic  correlation  effects  and  banding  effects  are  appar¬ 
ently  absent  in  our  earlier  one-electron  model. 


In  the  one-electron  model  one  can  approximate  U  by  a  differ¬ 
ence  between  two  terms,  the  direct  hole-hole  Coulomb  repulsion  F 
and  static  relaxation  term  R,  thus  U  =  F  -  R.  R  accounts  for  the 


discrepancies  are:  (l)  the  absence  of  the  high  energy  shoulder 
at  511  eV  in  the  0  KW  theoretical  lineshape,  and  (2)  the  sig¬ 
nificant  underestimation  of  the  low  energy  peak  at  50  eV  in  the 
Si  L23W  lineshape. 

That  good  agreement  between  the  cluster  LCAO-MO  theory  and 
experiment  should  be  obtained  is  perhaps  not  surprising.  The 
observed  Auger  lineshapes  in  the  more  covalent  Si02  are  very 
similar  to  the  corresponding  lineshapes  in  more  ionic  solids 
where  the  cluster  LCAO-MO  model  is  expected  to  be  valid.  Thus, 
for  example,  the  S  and  P  L23W  Auger  lineshapes  obtained  from  the 
highly  ionic  Li  salts  are  very  similar  in  all  aspects  to  the  Si 
L23W  lineshape  (3).,  Furthermore,  the  0  KW  lineshape  from  the 
more  ionic  MgO  is  essentially  identical  to  the  0  KW  lineshape 
in  Si02  (4).  The  valence  electron  energy  levels  for  Si02  have 
been  probed  by  photoemission,  by  Si  K.  and  L0  ,  x-ray  emission, 
and  by  0  Ka  x-ray  emission  (5).  These  data  have  been  well  des- 
cribed  within  the  Si04  cluster  one  electron  LCAO-MO  model  (l); 
indeed,  in  our  work,  the  orbital  energies  were  obtained  from  the 
XPS  data  and  the  LCAO-MO  populations  were  adjusted  to  reproduce 
the  XES  data.  Furthermore,  the  atomic  Auger  matrix  elements 
giving  the  Auger  intensity 
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state  holes  in  the  Auger  process)  were  not  evaluated  in  our  theory 
but  were  estimated  from  gas  phase  Auger  data  (2). 


Despite  these  successes  of  the  M.O.  model,  one  can  justifi¬ 
ably  question  the  validity  of  utilizing  the  cluster  one-electron 
M.O.  model  to  study  Auger  transitions  in  covalent  solids.  The 
Auger  process  produces  two  final  state  holes,  whereas  the  x-ray 
emission  and  photoemission  processes  produce  just  one  hole.  In 
general,  these  two  final  state  holes  see  one  another,  a  repulsive 
interaction  U  exist  between  them,  and  their  relative  motion  is 


correlated.  This  interaction  exhibits  itself  in  the  energy  of 
the  Auger  electron. 
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where  the  E's  are  the  one-electron  binding  energies  of  core  and 
valence  electrons,  respectively.  In  Si02,  U  is  substantial,  vary¬ 
ing  from  5  to  20  eV  depending  on  the  central  atom  and  the  localized 
nature  of  the  final  states  involved.  Hole-hole  correlation  may 
also  substantially  alter  the  Auger  intensities.  We  note  that  the 
hole-hole  correlation  effects  on  the  intensities  of  a  single  atom 
are  in  some  sense  already  included  in  our  model  because  the  experi¬ 
mental  atomic  Auger  matrix  elements  were  utilized.  Thus,  only 
the  interatomic  correlation  effects  and  banding  effects  are  appar¬ 
ently  absent  in  our  earlier  one-electron  model. 

In  the  one-electron  model  one  can  approximate  U  by  a  differ¬ 
ence  between  two  terms,  the  direct  hole-hole  Coulomb  repulsion  F 
and  static  relaxation  term  R,  thus  U  =  F  -  R.  R  accounts  for  the 


shift  in  the  binding  energy  of  one  hole  state  orbital  due  to  the 
presence  of  the  other  hole.  In  the  M.O.  model,  F  and  R  are 
expanded  in  the  atomic  equivalents 
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where  the  F°  are  the  atomic  Slater  integrals  between  atomic 
orbitals  i  and  j,  and  r^  is  an  atomic  intra-shell  relaxation 
energy,  both  defined  previously  (l). 

A  problem  with  the  cluster  one-electron  M.O.  model  is  that 
F  and  R  are  cluster  size-dependent.  Both  F  and  R  go  to  zero  as 
the  cluster  size  increases.  Thus  the  question  arises,  what  is 
the  appropriate  cluster  size?  in  the  one-electron  band  model  for 
a  solid,  F  and  R  are  of  course  zero.  This  is  merely  consistent 
with  the  notion  that  the  holes  are  completely  free  to  roam  about 
the  solid  or  cluster,  and  hence  are  completely  delocalized  from 
one  another.  There  is  experimental  evidence  (6-8)  that  in  some 
metallic  solids,  such  as  in  Cu,  Zn,  and  Ni,  the  final  state  holes, 
even  in  the  conduction  band,  are  not  free  to  roam  about  the  solid, 
but  are  forced  to  remain  on  the  atom  in  which  they  were  created 
(the  atom  with  the  initial  core  hole) .  This  is  evident  because 
the  experimental  Auger  intensity  does  not  reflect  the  density  of 
states  (DOS)  of  the  band  but  rather  exhibits  multiplet  structure 


consistent  with  an  atomic  model.  Furthermore,  the  apparent 
magnitude  of  U  is  consistent  with  the  hole-hole  interaction  one 
would  expect  if  both  holes  were  localized  on  the  same  atom. 

Electron  or  hole-hole  correlation  in  metallic  solids  has 
been  considered  within  the  context  of  the  Hubbard  (9)  or  Ander¬ 
son  ( 10)  model.  Within  these  models,  Cini  (ll)  and  Sawatsky  (12) 
have  independently  examined  the  role  of  electron  correlation  in 
metals,  and  have  attributed  the  behavior  indicated  above  to  hole- 
hole  correlation.  More  specifically,  if  the  effective  hole-hole 
interaction  U  is  large  compared  to  the  band  width  T,  energy  con¬ 
servation  forces  the  two  holes  to  remain  localized  on  the  atom 
where  they  were  created,  producing  an  Auger  intensity  consistent 
with  the  atomic  model. 

in  this  work  we  wish  to  consider  hole-hole  correlation  within 
the  cluster  LCAO-MO-CI  model  with  a  parameterized  Hamiltonian. 
Electron  correlation  in  molecules  is  normally  considered  within 
the  context  of  configuration  interaction  theory.  it  has  the  advan¬ 
tage  of  being  easily  applied  within  the  cluster  model  where  the 
two-center  Coulomb  repulsion  integrals  are  easily  included.  More 
importantly,  it  can  be  easily  applied  to  ionic  and  molecular 
insulators  and  more  covalently  bonded  systems  such  as  SiC>2,  which 
may  involve  more  than  one  element,  and  hence  involve  several  dif¬ 
ferent  interaction  energies.  In  Section  2  we  summarize  the  general 
theory,  providing  working  formulas  for  implementing  the  Cl  theory. 
To  elucidate  the  effects  of  hole-hole  correlation  and  attempt  to 


understand  why  the  one -electron  cluster  molecular  orbital  model 
may  be  adequate  for  interpreting  Auger  lineshapes  in  some  co¬ 
valently  bonded  solids,  we  examine  the  role  of  Cl  mixing  on  a 
model  problem. 

In  Section  3  we  present  numerical  results  appropriate  for 
the  0  KW  and  Si  L23VV  Auger  lineshape  in  Si02.  An  attempt  is 
made  to  account  for  the  two  described  discrepancies  between  the 
theoretical  and  experimental  Auger  lineshapes.  The  effects  of 
cluster  size  and  the  two-center  Coulomb  integrals  are  also 
examined . 


II.  CLUSTER  LCAO-MO-CI  THEORY  WITH  PARAMETERIZED  HAMILTONIAN 
1.  General  Formulations 

The  LCAO-MO-CI  theory  and  its  application  to  molecules  is 
well  known  (13).  We  briefly  summarize  here  the  general  formula¬ 
tions,  also  making  the  necessary  definitions  and  indicating  our 
approximations  and  limitations. 

The  molecular  orbitals  are  constructed  by  a  linear  combina¬ 
tion  of  atomic  orbitals. 
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in  the  usual  manner,  the  coefficients  determined  from  the  solution 
of  the  secular  equation. 
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The  atomic  orbitals  at  this  point  are  completely  unspecified; 
they  may  be  Slater-type  atomic  orbitals,  Hartree-Fock  atomic 
orbitals,  or  even  symmetry  adapted  atomic  orbitals  involving  more 
than  one  atomic  center.  They  are  specified  through  the  choice  of 
orbital  matrix  elements. 
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In  this  work,  the  matrix  elements  will  be  treated  as  parameters 

such  as  the  tight  binding  parameters  in  a  band  calculation  ( 14) ; 

that  is,  the  and  are  adjusted  so  that  the  resultant  one- 

2 

electron  M.O.  energy  levels  and  orbital  populations  c  ^  are 
similar  to  some  previously  determined  one-electron  density  of 
states  (DOS) .  We  have  reported  previously  our  best  estimates  of 
the  DOS  and  orbital  populations  in  Si02  (l).  They  were  determined 
from  experimental  data  and  ab-initio  cluster  M.O.  calculations. 

The  decision  here  to  set  the  off-diagonal  overlaps  equal  to 
zero  is  a  matter  of  convenience  for  this  model  problem.  In  general, 
they  are  not  zero;  however,  it  has  been  shown  (15)  that  by  includ¬ 
ing  them,  one  in  effect  scales  the  V. ..  Since  our  V. .  are  para- 
meterized  anyway,  the  ultimate  effect  of  ignoring  the  off-diagonal 
overlaps  is  minimal.  We  have  discussed  elsewhere  (16)  the  effect 
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of  on  the  distribution  of  local  versus  bonding  charge  and  the 
resultant  effects  on  the  Auger  intensity.  These  effects  are  not 
important  in  this  work. 

Within  the  configuration  interaction  theory,  we  shall  limit 
our  considerations  to  filled  bands.  This  means  that  all  configura¬ 
tions  describing  the  system  after  the  Auger  decay  will  have  exactly 
two  holes,  and  we  can  write  the  final  state  wavefunction  as  a  sum 
of  the  two  hole  configurations, 

^vv(fl'?2)  =  Eas  Dw',ap  <Pa(*l)  (8) 

Here,  the  <p^(r)  are  the  molecular  orbitals  spanning  the  cluster 
as  constructed  above.  The  restriction  to  filled  bands  is  a  good 
approximation  in  an  insulator  such  as  Si02,  as  the  unfilled  con¬ 
duction  bands  are  approximately  10  eV  above  the  filled  valence 
bands,  and  thus  an  insignificant  amount  of  configuration  mixing 
between  these  two  bands  should  occur  for  realistic  conditions. 

This  approximation  becomes  less  applicable  for  semiconductors  and 
metals  where  the  band  gap  is  small  or  zero. 

The  mixing  coefficients  ^  in  eq.  (8)  are  determined  by 

the  secular  equation 

)  ^w*  ys  ~  0  ^ 

where  the  full  Hamiltonian  now  includes  both  one  and  two  electron 
terms  H  =  h^  +  h2  +  U12.  Consistent  with  eq.  (7),  the  matrix 
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W 


elements  are 


H rj  =  <(  I  H  I  ^ ^  *  i~y 

'  < ^ ft ,  I l^xi  $s)  £or  +  U n 


(10) 


and 


ciCt>  J  ^  ^  ^  s  S 


(11) 


where 


=  1^12-1  ~  2-m» 


Ccir)-.  C(\/i  aXin'2n  ^  »*  ^  ■  (12) 


In  eq.  (12),  the  zero  differential  overlap  (ZDO)  approximation  has 

been  employed  to  reduce  the  two  center  integrals  U  Q  .  to  a  sum  of 

ap ,  y  6 

atomic  Coulomb  integrals  umn.  Throughout,  we  use  the  notation 
- 1  2—12 

<12|ri2  |  12>  =  [l  I r^2  |  2  ]  to  indicate  the  two  electron 
integrals. 

Those  integrals  involving  two  centers  are  evaluated  using  an 
approximation  due  to  Mataga  and  Nishimoto  (17) 

2 


u  =  e 
mn 


rR  . 
“  lun+V  J 


(13) 


and  the  one  center  integrals 


%  -  C*„2  I  r12  I  O  -  «„ 


(14) 
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are  treated  as  parameters.  In  eq.  (13),  Rmn  is  the  intemuclear 
distance  between  the  two  centers,  and  e  is  the  charge  on  the 
electron.  This  approximation,  found  in  many  semi -empirical  M.O. 


programs,  simply  takes  an  appropriate  average  of  two  extremums; 

2 


for  large  =  e  /R^,  for  very  small  R^,  =  (u^uJ/2, 

where  Un,um  are  the  one  center  integrals.  rR  in  eq.  (14)  is  the 
atomic  static  relaxation  term  as  mentioned  previously  (see  eq.  4); 
the  one  center  Coulomb  matrix  element  can  be  approximated  from  the 


tabulations  of  Mann  (18).  In  practice,  we  have  treated  u  as  a 

mn 


parameter  observing  hole-hole  correlation  effects  as  u  is 

mn 


increased. 

The  Auger  electron  energy  is  now  given  in  terms  of  the  eigen¬ 
values  of  eq.  (9) 


Ecw'  Ec  “  Ew ' 


(15 


in  contrast  with  eq.  (2)  which  involved  the  one-electron  M.O. 
energies.  The  Auger  intensity  now  involves  the  two  hole  DOS 


lew' s  T  K-fcWi  ni'l  Yvv.  yi 


vv 

2 


^  1 2  Dvv>ft  c*a  cfcbl  l^cab 


(16 


where  we  have  inserted  eqs.  (l),  (5),  and  (8),  and  made  the  usual 
assumption  that  only  the  one  center  Auger  matrix  elements  are 
significant  (l).  The  indices  "a"  and  "b"  label  atomic  orbitals 
on  the  atom  with  the  initial  core  hole;  thus  the  Auger  lineshape 


reflects  the  two  hole  local  DOS.  It  is  this  feature  which  hastens 
the  convergence  of  the  cluster  expansion  result  to  the  bulk  limit 
and  thus  provides  a  reasonable  probability  for  the  success  of  the 
cluster  expansion  in  this  work. 

The  above  analysis  ignores  the  role  of  hole  (electron)  spin. 
Actually,  the  spin  symmetry  decreases  the  number  of  two  hole  con¬ 
figurations  which  have  to  be  included  in  the  final  state  wave- 

function.  Configurations  of  the  "diagonal"  type  tp^tp^  only  take 

.1# 

the  singlet  form  2e(ap-3a).  Non-diagonal  configurations  can 

take  both  the  singlet  and  the  triplet  forms 


-V;  K ue-e.) 


2K  (K$s'  ^ 


e>a 


singlet 


triplet 


Since  singlet  final  states  can  couple  only  to  singlet  initial 
states  (and  similarly  for  triplet  states) ,  eq.  (16)  has  the  form 


J/s 


r  ,  r  ,  o  ( i  for  the  diagonal  terms,  and 

{"'cab  / 


Dv„'j(*b  [  cSb  T  caaJ(Mcal,*  Mcbn)  (17) 

for  the  off-diagonal  terms,  where  the  +/-  sign  denotes  the 
singlet/triplet  states,  respectively.  For  a=b,  the  diagonal  and 
singlet  off  diagonal  terms  are  the  same  except  for  the  additional 
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\[2  factor  in  the  singlet  off-diagonal  term,  the  triplet  terms 
are  zero.  For  a^b,  the  triplet  terms  are  generally  smaller  than 
the  singlet  terms.  The  spin  symmetry  has  the  effect  of  reducing 
the  size  of  the  secular  determinant  since  the  singlet  and  triplet 
configurations  form  non-interacting  sub- determinants. 

A  common  problem  with  the  Cl  technique  is  the  rapid  increase 
in  the  number  of  configurations  required  with  increasing  molecule 
or  cluster  size.  The  local  nature  of  the  two  hole  DOS  reflected 
in  the  Auger  lineshape,  the  relative  unimportance  of  the  unfilled 
bands  in  insulators,  and  the  spin  symmetry,  all  points  made  above, 
assist  us  greatly  in  making  the  problem  tractable,  in  this  work 
the  number  of  configurations  never  exceeded  50;  however,  this  was 
in  some  cases  the  minimum  number  necessary  for  achieving  our  semi- 
quantitative  results. 

2.  Model  Problem. 

To  elucidate  the  role  of  configuration  mixing  and  hole-hole 
correlation  in  Auger  lineshapes,  we  examine  here  and  solve 
analytically  a  four-orbital  problem.  As  indicated  above,  the 
limitation  to  four  orbitals  and  four  or  less  configurations  is 
a  severe  restriction,  but  it  does  retain  the  important  features 
and  yet  is  simple  enough  to  offer  an  intuitive  feeling  for  the 
correlation  effects.  Indeed,  this  is  an  important  aspect  the 
cluster  LCAO-MO-CI  model  has  over  the  Hubbard  model  which  treats 
infinite  systems  (11,12). 
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The  restriction  to  a  four-orbital  problem  and  the  desire 
to  have  at  least  two  identical  sub-clusters  limits  us  in  this 
model  problem  to  a  linear  system  with  triatomic  sub-clusters 
B-A-B,  where  we  assume  that  the  orbitals  on  B  are  hybridized 
(e.g.,  sp  hybridized).  In  this  respect,  B  will  play  the  role 
of  the  Si  and  A  of  the  0  atoms  in  Si02,  as  it  is  generally  agreed 
(19, 20)  that  Si  is  highly  sp  hybridized,  whereas  0  undergoes 
little  hybridization.  Table  I  summarizes  two  cases  which  are 
within  these  limitations,  case  a  involving  bonding  and  anti¬ 
bonding  orbitals  on  A  and  B,  case  b  involving  nonbonding  orbitals 
on  atom  A.  The  secular  determinant  is  most  general  for  case  a; 
case  b  being  the  special  case  when  aQ  =  and  V  =  T. 

2.  (Case  a) 

Solution  of  eq.  (6)  for  case  a  gives  the  four  M.O.'s 


I  i/li)  -  (  #A,+  a~  t  h,:  ~  Q  ) 

hi, 7^.  A,“  tin  ) 


(18) 
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and  the  M.O.  energies  are 

6W»i  ~ 

l  +  Vau 

(20) 

£  lii/iv  - 

F(F 

r+yi  • 

Appropriate  for  the  case,  r«V,  eq.  (20)  has  been  written  to 


emphasize  the  existence  of  two  closely  spaced  bonding  orbitals 
eb-,  and  two  closely  spaced  antibonding  orbitals  ea±.  In  an 
extended  system  one  would  obtain  a  "band"  of  energies  at  and 
one  at  ea#  with  band  widths  of  the  order  of  T  and  a  band  gap  of 
the  order  V. 

Assuming  in  "a"  that  the  bonding  orbitals  are  filled,  and 
the  antibonding  orbitals  empty  and  sufficiently  removed  from  the 
bonding  orbitals  (v  >  u) ,  we  have  just  three  possible  two-hole 
configurations:  liijti)}  |  i;  ii  ^  .  Further  simplifi¬ 

cation  results  from  symmetry  requirements  which  exclude  any 
mixing  of  theli,ii>  configuration  with  the  two  diagonal  configura¬ 
tions.  The  problem  reduces  to  a  2  x  2  problem  with  the  secular 
determinant 


<v  I 


•  v) 

<?£*  * 


''I'1) 


a  ~T~ 


z  z 


L 


+  - 


Zej+  l|,  (All)"  4 


(21) 
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where 


(22) 


and 


MAH)1-*  [(V^w'U 

=  L|Nt  [uA,A1*(o:tQ?;)l«nini+ 


v) 

aV(M^  +  MM()  ] 


(23) 


represent  respectively  the  effective  interaction  between  two  holes 
on  a  single  B-A-B  cluster  and  on  two  different  clusters  (b1~A^-b2) 
and  (b2-A2-B3)  in  the  four  atom  system.  Solutions  of  the  secular 
equation  give  the  eigenvalues 


e ■  +u,(Mir+u,(AHr 

+  uIi(ah)'h>ui1(ah',"']  h 
i  Ml(AHr-  U,(AWft 

-uu(M)"Th  ♦fu,(Mr-«i»i*»r],J 
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_ 


■lil 


and  mixing  coefficients 

-D,j  .  2Er-f^-(U,(M<r-^^a(AHr) 
■D'vi'  U,(aH^-.uu(AvO+' 


The  energy  of  the  |i,ii>  configuration  can  be  given  by  the 
expression 


E„  =  <  2'4  (  ./»  +  » ,  i  )  |  fi  I  (  i »  *  <>j  <  )  ) 


-  ^  +  M,  ( AH)* 


(26) 


It  has  already  been  assumed  that  U  «  V.  To  reduce  the  complexity 
of  eqs.  (24)  through  (26),  we  shall  assume  here  that  T  «  U^AH). 
With  this  assumption,  the  first  term  in  brackets  inside  the  square 
root  in  eq.  (24)  is  negligible  compared  to  the  second.  Furthermore, 
we  define  the  average  quantities  eb  =  l/2(eb  +  eb) ,  U-^AB)  = 
1/2(U1(AB)++  +  U1(ab)__),  etc.  Then  eqs.  (24-26)  become 

e"  =  2  €b  +  u12(ab) 

E°  “  2  eb  +  ui^AB) 

E+  =  2  eb  +  U1(AB)  (27) 


and 
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5Eu3H=5=3 


r  » 


(£)*  '  ■  - 


Inserting  eq.  (28)  into  eqs.  (8)  and  (18)  gives 

^  -  Y?  N[($Ai  +  q  ^HzV2) 

+  (^AZ+q 


(28) 


(29a) 


7  ( ‘-Paz  +' ^ 4^ng')(i's  ($a2  +  a$tt2)(2) 


(29b) 


where  we  have  defined  Na  =  l/2(N+a++N  a  ).  It  is  clear  that  in  the 
+  ° 

states  ^  and  if  the  two  holes  are  localized  on  the  same  B-A-B 
cluster;  in  the  ii  state,  they  are  localized  on  different  clusters. 

If  we  insert  eq.  (29)  into  eq.  ( 16) ,  we  can  obtain  either  the 
A  or  B  Auger  intensity.  Assuming  the  initial  core  hole  was  on 
atom  A^,  the  Auger  process  picks  out  the  DOS  local  on  atom  A^,  or 
<Pai<Pai  terms*  If  the  initial  core  hole  was  on  atom  B2#  it  picks 
out  the  DOS  local  on  atom  B2,  or  the  terms  ^1^112 '  ani^ 

^H2^H2*  The  results  are 
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^A|  “  ° 


-^A|  "  I  ^cAA  I  "2  -^A  (30) 

*i\  -  ^^iMCAAr=  t 


-  i N'° 1 

1  Mess  "  Mcpfa  ( 

4=  2nV 

1  Mc5  p  ) 

(31) 

1 82  =  iNV 

1  M epp  1 

where  IA  is  the  total  intensity  for  the  "band"  and  M  is  the 
atomic  Auger  matrix  element  (eq.  (l))  involving  the  core  orbitals 
on  atoms  A  and  B  as  appropriate.  Although  we  have  not  specifically 
introduced  s  and  p  orbitals  on  atom  B  to  this  point,  we  assume  here 
for  convenience  that  the  orbitals  h^  and  h2  are  sp  hybrids,  so 
that  s(p)  =  yy-  (h^t  h2) . 

The  intensities  in  eqs.  (30)  and  (31)  can  be  contrasted  with 
those  obtained  from  the  fold  of  the  one  electron  DOS.  They  are 

it*  N-’|McAA|%  I1 

^ At  -  lwCMr  *  i rA  02) 

* *,  5  n*-*  1  m£AA  r « 
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I.-  H 


v< 

1  az  ~ 

(q  )  |  M  Cpp  j 

-r-  1  ^  '' 

,  •  Z  |  j  z 

-BZ  = 

2  N.N+  (q  cT)  |  | 

r 1,1  = 

B2 

w.2  r 

energies 

E“jli  - 

2s.k-f  (U,(A»r+  L4,X(AH)_ 

E  '' ' '  i 

*  U "  C  A  H  i 

(33) 


(34) 


E  1/1  =  2<  *  (U,(W'f  t  Ulz(A«r+)/2 


The  intensities  in  eqs.  (3o)  and  (3l)  are  appropriate  when 
r« l^(AB)«V;  those  in  eqs.  (32)  and  (33)  when  U1(ab)«T«V. 

A  comparison  of  these  results  reveals  some  interesting 
features.  For  the  case  of  the  initial  core  hole  on  atom  A,  eqs. 
(27),  ( 30)  and  (32)  reveal  an  overall  shifting  of  the  same  total 
intensity  from  the  term  involving  non-local  final  state  holes 
(one  hole  on  B^-A^-I^;  the  other  on  B2_A2~B3^  to  the  terms 
involving  local  final  state  holes  (both  holes  on  B2~A^~B2)  as 
U  increases.  This  is  a  general  result,  one  which  will  be  examined 
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numerically  for  larger  clusters  in  Si02  in  Sec.  3.  The  situation 
is  significantly  different  when  the  initial  hole  is  on  atom  B. 
First,  as  increases,  the  total  intensity  changes  due  to  increas¬ 
ing  interference  effects  between  the  M  and  M  contributions 

CSS  CPP 

in  the  large  U  case.  Second,  even  when  U  »  T,  the  intensity  from 
the  term  involving  non-local  final  states  holes  i“  does  not 
necessarily  go  to  zero.  In  fact,  it  may  be  larger  than  I*  •  this 
depends  on  the  signs  of  Mcsg  and  MCpp*  In  summary,  we  can  say 
for  U  »  r,  the  Auger  process  on  the  A^  atom  samples  the  DOS  local 
to  the  smaller  sub-cluster;  the  Auger  process  on  the  B2 

atom  samples  the  DOS  on  the  larger  Bi-A1-B2_A2"B3  cluster  but  with 
significant  hole-hole  "correlation"  effects. 

An  examination  of  several  cases  reveals  that  quite  generally 
I+B  is  less  than  I  B  ;  i.e.,  MCpp  and  Mcgg  have  opposite  signs. 

For  the  case  at  hand,  the  sp  hybrid  and  a  Is  initial  core  hole 

z 

on  B2 #  we  have 

(N'css  J  W'cpp)1--  |  <ls  + 

00-7  l  (35) 

-  \  R(s^ppV)  --i  RVsJpp) 

where  R  is  the  radial  integral 


W  ^is  Rut  Jji  r r,2  • 

O  &  •  Li  I 
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For  atoms  in  the  first  rows  of  the  periodic  table,  R(ssss)  is 

O 

positive,  R(sspp)  is  negative  (21).  For  the  K-LL  transitions 
in  atoms,  Asaad  (22)  has  found  that  Cl  effects  decrease  the  KL^L^ 
intensity  and  increase  the  KL^L^  intensity.  McGuire  (23) 
determined  that  Cl  effects  decrease  the  L23M1M1  intensity  in  Cu 
and  Zn  and  increase  the  L23M23M23  intensity. 

2.2  (Case  b) 

Case  b  involves  four  nonbonding  orbitals  on  four  A  atoms  as 

indicated  in  Table  I.  It  is  a  special  case  of  "a",  namely  when 

e  =e  and  V  =  T;  hence  the  molecular  orbitals  and  energies  are 
O  H 

still  given  by  eqs.  (18)  to  (2o).  More  specifically,  these 
equations  give 

|i/i,)  =  q  ~  ) 

(3  7 

|  ill  /  iv)  ~  hl-r  (  \\  ~  Q  "  ^ A2  “  3  ^ 


?  ^  J-  4-  {I. 

a  -  -  l 


a  + 


A  2 


(f5  ±  \) 


e  ‘in/ iv  '  ^  A  4  ^  (-is  t  i) 
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Assuming  that  all  the  nonbonding  orbitals  are  filled,  we 
have  ten  possible  two-hole  configurations.  In  the  event  U  «  T, 
a  negligible  Cl  mixing  occurs,  and  the  two  electron  energies  and 


where  the  initial  core  hole  is  assumed  to  be  on  atom  A2  and  the 
proper  choice  in  eq.  (41)  is  obvious.  The  total  sum  E  I  equals 


In  the  event  U  T,  the  analytical  solution  of  a  10x10  secular 
determinant  is  formable  even  though  the  problem  is  simplified 
somewhat  by  symmetry.  Much  more  instructive  at  this  point  is  to 
examine  the  problem  when  U  »  I\  This  can  be  examined  more  simply 
by  omitting  the  construction  of  the  M.O. 's  and  constructing  the 
two-hole  configurations  directly  in  terms  of  the  atomic  orbitals 
<PA«  The  Hamiltonian  matrix  becomes 
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where  the  |ij>  now  refer  to  atomic  configurations  <pa  .<p  . .  in 

**  1  nj 


eq.  (42)  we  have  omitted  the  2  term  from  each  of  the  diagonal 
elements  and  utilized  the  ZDO  approximation  as  usual.  The 


symbols  u—  refer  to  Coulomb  matrix  elements  on  the  same  center 


and  between  the  1st,  2nd  and  3rd  nearest  neighbors,  respectively. 
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The  diagonalization  of  eq.  (42)  produces  the  same  eigen¬ 
values  and  eigenvectors  as  the  diagonalization  of  the  Hamiltonian 
matrix  constructed  from  the  ten  M.O.  configurations.  In  general, 
the  Hamiltonian  matrix  is  more  difficult  to  diagonalize  in  the 
atomic  basis  because  symmetry  considerations  do  not  simplify  the 
problem  as  in  the  M.O.  basis.  Nevertheless,  if  u^>>r  and  thus 
also  j  )  >:>  r/  this  matrix  is  essentially  already  diagonal 

and  the  solution  becomes  trivial.  The  Auger  process  samples  the 
DOS  local  on  atom  A2;  the  entire  Auger  intensity  |mcAA|2  coalesces 
to  one  narrow  peak  at  energy  2a  +  u.  . 

The  case"b"  as  described  is  appropriate  for  the  non-bonding 
orbitals  on  0  in  the  Si02  system.  In  Sec.  3  we  will  examine  the 
coalescence  to  a  single  peak  as  increases  relative  to  r,  and 
also  examine  the  effects  of  cluster  size  and  the  effects  of  the 

two  center  u.  . . 

ID 

III.  RESULTS  AND  DISCUSSION 

We  have  performed  various  M.O.  calculations  on  several  differ¬ 
ent  clusters,  always  with  a  minimal  basis  set  and  at  times  with  a 
partial  minimal  set.  The  valence  one  electron  DOS  in  Si02  fall 
roughly  into  three  major  group J:  the  "02g"  bonding  states  grouped 
about  20  eV  below  che  edge  (the  4a.^  and  3t2  orbitals  of  the 
tetrahedral  Si04  cluster) ,  the  "02p"  (b)  bonding  states  grouped 
about  6  eV  below  the  edge  (the  5a^  and  4t2  orbitals) ,  and  the 
"°2p"  nonbonding  states  located  just  below  the  edge  (the  le. 
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5t 2  and  lt^  orbitals)  .  We  have  examined  the  02p(b)  and  tlie 
02p(nb)  states  on  the  Si403  and  Si807  clusters  with  0  as  central 
atom,  the  0 2g/  °2p(k)  and  02p(nb)  states  on  a  Si504  cluster  with 
Si  as  central  atom,  and  the  02p(nb)  density  of  states  on  03,  07, 
and  clusters.  A  beta  crystobalite  crystal  structure  (Si-O-Si 

angle  equal  180°)  is  used  throughout  with  Si-0  bond  length  equal 

O 

1.5  A  (24).  The  initial  core  hole  is  placed  in  a  central  Si20 
sub-cluster,  nearest  neighbor  sub-clusters  are  placed  in  the 
proper  crystal  structure  building  outwardly  to  assemble  the  com¬ 
plete  cluster. 

The  matrix  elements  in  eq.  (6)  were  treated  as  parameters, 
determined  so  that  the  resultant  one  electron  M.O.  energies  and 
orbital  populations  are  similar  to  our  previously  determined  Si02 
DOS  (l).  An  optimal  set  of  matrix  element  parameters,  determined 
from  the  above-indicated  calculations,  is  reported  in  Table  II 
and  compared  with  corresponding  tight  binding  parameters  as 
recently  reported  (14,25).  In  general,  the  agreement  is  quite 
satisfactory  considering  our  calculations  involve  a  minimal  basis 
set  on  rather  small  clusters.  The  v(°2S'°2s^  an<^  V^°2p'°2p^  e^-e“ 
ments  include  both  direct  and  indirect  "hopping"  or  interaction 
terms,  i.e.,  hopping  terms  which  involve  the  Si  atom  are  included. 
In  fact,  the  presence  of  the  Si  atoms  in  the  0^  02p(nb)  calcula¬ 
tion  is  exhibited  only  through  these  terms.  Only  hopping  matrix 
elements  involving  two  oxygen  atoms  connected  through  a  common  Si 
atom  are  assumed  to  be  significant,  all  others  are  set  to  zero. 
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Likewise,  only  nearest  neighbor  bonded  Si-0  elements  are  included. 

Appropriate  values  for  the  matrix  elements  involving  the  Si  sp3 
hybrid  orbital  h^  are  obtained  as  indicated  in  the  table. 

The  values  of  u  tabulated  in  Table  II  are  obtained  from  the 
integrals  of  Mann  (18);  however,  a  static  relation  term  r  has 
been  subtracted  off.  These  have  been  estimated  previously  (l)  to 
be  6  eV  for  02s  and  2p  orbitals,  and  2(1)  eV  for  Si  3s (3p)  orbitals 
Huang  (15)  has  given  the  formula  for  obtaining  the  u(Sih)  result 
from  the  -u(Si3s/p)  elements.  u(o2p)  has  a^-so  been  determined 
empirically  by  Oleari  et.  al.  (26)  who  obtained  a  value  of  14.5  eV. 

The  empirical  result  is  said  to  include  correlation  effects  (i.e., 
static  relaxation  effects) . 


1.  0  KLL  Line shape 

Results  for  the  0  KLL  Auger  lineshape  obtained  from  the 
linear  Si403  cluster  involving  the  "02p"  bonding  and  nonbonding 
orbitals  are  given  in  Fig.  2.  These  results  were  obtained  from  a 
12-orbital  system;  namely  3  Si^  =  l/\[2  (h^-h2)  orbitals,  3  02p(b)  , 
and  6  02p(nb)  atomic  orbitals.  In  this  system,  symmetry  allows 
mixing  only  between  the  02pz  and  Si^  orbitals  giving  M.O.'s 
belonging  to  the  02p  bonding  and  antibonding  groups;  however, 
the  antibonding  M.O.’s  are  unoccupied.  Thus  the  b-b,  b-nb,  and 
nb-nb  two-hole  configurations  are  mixed  in  separate  Cl  calculations 
The  results  show  clearly  the  effects  of  increasing  correlation 
(increasing  u) ,  most  notably  in  the  coalescence  of  most  of  the 
Auger  intensity  into  a  rather  narrow  band  around  502  eV.  Clearly, 
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this  band  corresponds  to  both  holes  localized  on  to  the  same 
Si-O-Si  cluster.  The  02p(nb)  lines  are  shifted  down  by  ~15  eV 
as  they  are  localized  on  a  single  0  atom;  the  02p(b)  lines  just 
11  eV  (the  approximate  hole-hole  repulsion  on  the  Si-O-Si  cluster 
in  the  09  (b)  M.O.).  Note  the  appearance  of  residual  non-local 
contributions  shifted  down  by  just  4  eV  (the  approximate  hole-hole 
repulsion  when  the  holes  are  on  different  Si-O-Si  clusters). 
Ignoring  these  residual  contributions,  the  Auger  02p  lineshape 
narrows  from  ~9  eV  to  ~2  eV,  giving  essentially  an  atomic  0 
Auger  spectrum. 

Since  the  greatest  share  of  the  502  eV  peak  of  the  0  KLL 
Auger  lineshape  comes  from  the  (nb,nb)  transitions,  we  can  for 
the  moment  ignore  all  other  contributions.  This  allows  for  a 
study  of  much  larger  clusters  and  hence  also  a  study  of  cluster 
size  effects.  Results  of  calculations  on  the  01&  clusters  are 
given  in  Fig.  3,  where  just  the  Auger  (nb,nb)  contributions  are 
examined.  The  one  electron  DOS  and  the  two  hole  DOS  resulting 
from  a  Cl  involving  37  configurations  are  shown.  Similar  trends 
are  revealed,  namely  a  coalescence  to  a  local  peak  shifted  down 
in  energy  by  an  amount  -u^  and  several  residual  non-local  peaks. 

Figure  5  summarizes  results  of  calculations  on  O^,  0^,  and 
O,,  clusters  when  all  two-center  hole-hole  repulsion  integrals 
are  set  to  zero.  Plotted  in  Fig.  4  is  the  intensity  in  the 
residual  non-local  peaks  relative  to  the  total  intensity  (always 
normalized  to  100  in  these  calculations)  and  the  difference  in 
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energy  between  the  centriods  of  the  local  and  non-local  contribu¬ 
tions  vs.  u/T.  u/T  is  the  one-center  hole-hole  repulsion  on  the 


oxygen  atom  divided  by  the  one  electron  02p(nb)  band  width.  The 
results  plotted  depend  only  on  u/r  and  not  on  u  and  T,  individually, 
and  no  other  parameters  are  present  in  these  calculations.  As  one 
might  expect,  the  results  from  the  different  cluster  sizes  converge 
with  increasing  u/T;  for  large  u/T  the  Auger  lineshape  is  increas¬ 
ingly  dominated  by  the  local  peak  where  both  holes  are  on  a  single 
0  atom.  in  the  physical  range  of  interest,  u/T  «  2  to  4,  it  would 
appear  that  the  07  cluster  size  is  already  sufficient. 

Also  plotted  in  Fig.  4  is  result  of  a  calculation  on  the  0.c 

±b 

cluster  when  the  two-center  integrals  u^j  are  included  as  esti¬ 
mated  from  the  Mataga-Mishimoto  approximation,  eq.  (13).  in  these 
calculations,  the  u^  increase  with  increasing  u/T.  The  inclusion 
of  the  u^j  appears  to  decrease  the  effective  size  of  the  one-center 
u^,  as  one  might  expect.  However,  it  should  be  pointed  out  that 
the  full  effects  of  cannot  be  completely  compensated  for  by 

decreasing  the  effective  u^.  The  shift  down  in  absolute  energy  of 
the  total  Auger  lineshape  increases  when  the  u^j  are  included, 
but  a  decrease  in  u^  would  decrease  this  shift. 

The  difference  peak  at  511  eV  in  the  0  KLL  Auger  lineshape 
(Fig.  l)  is  estimated  to  have  about  20$  of  the  intensity  of  the 
theoretical  Auger  peak  at  504  eV  (the  shake  satellite  intensity 
is  not  included  in  the  504  eV  peak  here) .  If  this  difference  peak 
is  attributed  to  the  non-local  Auger  contributions,  the  7  eV 
experimental  separation  between  the  local  and  non-local  peaks 
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provides  an  estimate  of  2.5  for  u/r  (from  Fig.  4,  using  the  0.c 
nonzero  plot)  .  Furthermore,  the  17$  experimental  intensity 

ratio  between  the  difference  peak  at  511  eV  and  the  total  (504  eV 
peak  plus  511  eV  peak)  indicates  a  u/T  value  of  2.5  (also  from 
Fig.  4,  non-zero  plot).  That  the  two  estimates  are  in  agree¬ 
ment  is  most  gratifying  and  supports  rather  strongly  the  assumption 
that  the  511  eV  shoulder  is  indeed  due  to  the  non-local  contributions. 

Note  that  the  0, results  obtained  for  zero  u. .  would  give  con¬ 
ic  13 

flicting  values  for  u/T. 

The  value  of  2.5  for  u/T  is  a  reasonable  one.  The  best 
estimate  of  u(°2p^'  as  ?iven  in  Table  II,  is  15  eV.  This  indicates 
a  value  of  6  eV  for  T.  Band  calculations  (14,19,20,25,27,28)  have 
given  results  varying  from  1  to  5  eV  for  the  bandwidth  of  the 
02p  nonbonding  band  with  perhaps  the  most  reliable  theoretical 
result  around  5  eV  (20,27).  X-ray  emission  data  and  UPS  data  in 
Si02  suggest  a  value  around  4  eV  (27,29,30,31). 


We  can  suggest  two  possible  causes  for  our  1-2  eV  overestimate 
of  r.  The  determination  of  a  bandwidth  from  a  relatively  small 
cluster  calculation  is  nebulous  at  best.  The  O^g  cluster  calcula¬ 
tions  provide  a  spectrum  of  eigenvalues  the  breadth  of  which  pro¬ 
vides  our  estimate  of  the  bandwidth.  This  is  probably  the  best 
procedure  since  in  band  calculations,  T  represents  the  full  width 
of  the  DOS.  Secondly,  our  results  depend  on  the  values  of  the 

u...  This  is  indicated  in  Fig.  5  for  the  two  0. c  cluster  calcula- 
13  lb 

tion,  one  including  the  vu  ^ ,  the  other  not.  If  the  Mataga-Nishimoto 
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approximation  underestimates  the  u„  (we  will  present  some  evi¬ 
dence  later  that  this  may  be  the  case) ,  our  estimate  of  T  will 
be  too  large.  Evidently  as  ^  increases,  both  the  relative 
intensity  and  energy  difference  curves  move  to  higher  -u/T  and 
hence  indicate  a  smaller  T. 

It  is  of  interest  from  an  intuitive  point  of  view  to  note  the 
relative  location  of  the  two  holes  formed  in  the  Auger  process. 
Table  ill  contains  the  probabilities  P^j  for  hole  occupancy  on 
the  i^n  and  j*"*1  0  atoms  in  the  C>2p(nb)  band  of  Si(>2.  For  simpli¬ 
city,  these  results  come  from  an  07  cluster  calculation  rather 
than  the  cluster.  The  two  "local"  Auger  peaks,  peaks  4  and  5 
(peak  numbers  are  indicated  in  Fig.  3),  clearly  correspond  to  a 
high  hole  occupancy  on  the  same  0  atoms  as  anticipated  for  large 
u.  The  three  "non-local"  peaks  in  order  of  decreasing  energy 
can  be  described  as  follows:  peak  (3)  the  hole  occupancy  is 
predominantly  on  2nd  and  3rd  near  neighbors  with  a  significant 
probability  (~. 11)  for  both  holes  on  the  same  center;  peaks  (2) 
and  (l),  the  hole  occupancy  is  predominantly  on  nearest  neighbors 
but  with  significant  probability  for  occupancy  on  2nd  and  3rd 
near  neighbors  also.  The  Auger  intensity  is  directly  proportional 
to  P^1  (both  holes  occupying  the  atom  with  the  initial  core  hole). 
It  would  appear  the  Auger  energy  is  primarily  determined  by  the 
total  of  the  one  center  occupancy  terms  as  the  "non-local"  peaks 
are  ordered  in  energy  accordingly.  The  relative  occupancy  of 
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nearest  neighbors  and  2nd  and  3rd  near  neighbors  is  unimportant 
in  determining  energy  apparently  because  the  1/R  behavior  reduces 
u12  already  to  ~4  eV  (when  u1  is  15  eV) ,  and  the  higher  terms  (e.g., 
U13'  u25*  are  onlY  slightly  smaller. 

2.  Si  L23W  Line  shape 

Results  for  the  Si  L^W  Auger  lineshape  obtained  from  the 

tetrahedral  Si504  cluster  is  given  in  Fig.  6.  These  results  were 

obtained  from  a  16  orbital  system;  namely  four  Si,  orbitals  on 

h 

the  central  Si,  one  Sih  orbital  on  each  of  the  edge  Si  atoms,  and 
an  0  2s  and  2p  orbital  on  each  of  the  four  0  atoms.  The  0  2p 
and  2py  nonbonding  orbitals  were  ignored  in  this  calculation.  The 
resultant  one  electron  M.O.’s  are  named  in  the  T^  point  group 
4ai3s),  3t2(3p),  5a-j(3s),  and  4t2(3p),  where  the  dominant  contribu¬ 
tion  of  the  Si  atom  in  each  M.O.  is  indicated  in  parentheses. 
Although  the  M.O.  calculations  are  performed  with  sp  hybrids  (h^) 
on  the  Si  atoms,  solution  of  the  secular  equation  automatically 
provides  for  the  proper  linear  combination  of  these  hybrids  to 
give  back  either  an  s  or  p  orbital  as  dictated  by  symmetry.  Upon 
performing  the  Cl  with  increasing  u^,  the  orbitals  rehybridize. 

The  number  of  two-hole  configurations  necessary  for  an 
appropriate  hole-hole  correlation  treatment  can  again  be  kept 
to  a  minimum.  The  energy  spacing  between  the  4a^,  3t2  orbitals 
(the  02g  band)  and  the  5a^,  4t2  orbitals  (the  02p(b)  band)  is 
large  compared  to  the  effective  hole-hole  repulsion  in  the  Si504 
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cluster,  so  the  configurations  from  the  02s°2s'  °2s°2p  and 
02p-02p  bands  can  again  be  mixed  in  separate  Cl  calculations. 
Group  theoretic  arguments  can  be  used  to  simplify  the  problem 


still  further.  Thus  in  the  "diagonal"  cases  (either  °2s°2s  or 


°2p-°2p)'  we 

have 

the  configurations  leading  to  the 

T. 
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1  T  5  T 

Since  only  terms  of  the  same  symmetry  and  multiplicity  mix,  the 
complete  Cl  is  easily  tractable  (e.g.,  in  the  diagonal  cases, 
mixing  occurs  only  between  the  two  ^A  terms  and  only  between 
the  two  ^T2  terms) . 

Expressions  for  determining  the  Auger  intensity  are  given 
in  Table  III.  These  expressions  involve  the  Cl  mixing  coefficients 
D  obtained  from  eq.  (9),  the  atomic  orbital  coefficients  c  obtained 
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from  solution  of  eq.  (6),  and  the  atomic  Auger  matrix  elements  in 

L-S  coupling  •  (  2S+^l)  •  Expressions  for  evaluating  Mc£/( 

have  been  reported  by  McGuire  (23)  and  others  (22,32,33).  The 

Mcpp ^  3 s ^  #  MCpp^ '  Mcpp^lp^'  and  Mcss^S^  matrix  elements  are 

zero  in  this  work,  as  shown  in  eq.  (17)  and  sec.  II. 1.  The  five 

terms  in  spherical  symmetry  split  into  the  and  ^2  terms 

in  tetrahedral  symmetry.  The  intensity  splits  up  equally,  one- 

fifth  into  each  of  the  terms  regardless  of  symmetry  (34).  As 

determined  by  McGuire,  M  (^S)  and  M  (^S)  are  of  different 

css  cpp 

2  2 

sign  so  that  the  a^  +  t2  state  has  a  smaller  intensity  than  the 
2  2 

a^  -  t2  state. 

Figure  5  compares  results  with  no  hole-hole  repulsion  effects, 

with  hole-hole  repulsion  effects  but  without  Cl  mixing,  and  with 

complete  Cl  mixing.  The  turn  on  of  hole-hole  repulsion  without 

Cl  mixing  causes  additional  lines  from  multiplet  splitting.  The 

turn  on  of  Cl  mixing  then  coalesce  these  terms  into  "local"  and 

"non-local"  contributions  for  each  band  (0o-0-  ,  0_-0»  ,  and 

2s  2s  2s  2p 

°2p°2p)*  In  t^1^s  instance,  the  local  contributions  correspond  to 
both  holes  localized  on  an  Si20  sub-cluster  with  «  11  eV, 

and  the  non-local  contributions  to  the  holes  localized  on  differ¬ 
ent  Si20  sub-clusters  with  w  4  eV.  Because  the  initial  core 

hole  was  on  the  central  Si  atom  with  a  local  Si  population  in  all 
four  Si20  sub-clusters,  the  Auger  intensity  now  does  not  coalesce 
into  just  the  local  peak  as  in  the  0  KLL  case,  but  a  major  part  of 
the  Auger  intensity  remains  in  the  non-local  contributions. 
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nevertheless,  the  effects  of  Cl  mixing  appear  to  be  evident 


in  the  Si  I^VV  lineshape.  The  local  contribution  of  the  ^s^p 
band  at  ~54  eV  now  matches  up  more  with  the  50  eV  peak  in  the 
experimental  Auger  lineshape  of  Fig.  1.  The  .4  relative  intensity 
of  this  peak  also  is  of  the  proper  magnitude  to  account  for  the 
intensity  under  the  experimental  peak  since  the  theoretical  peak 
was  way  too  small  in  the  earlier  work.  The  shoulder  at  68  eV  in 
the  experimental  lineshape  we  have  previously  suggested  may  be 
due  to  a  shake  satellite.  It  now  appears  it  could  be  due  to  the 
local  contribution  of  the  °2p°2p  band  at  ~70  eV  with  a  relative 
intensity  of  .7.  The  intensity  of  .7,  somewhat  less  than  the 
nearby  1.2  intensity  of  the  contributions  at  63  eV,  seems  qualita¬ 
tively  correct  if  it  is  to  account  for  the  68  eV  shoulder.  The 
small  local  contribution  of  the  0~  0,,  band  at  37  eV  falls  under 
the  Si  L23  V  peak  which  we  examined  earlier  (l). 

It  appears  that  the  local  contribution  is  sufficiently  removed 
in  energy  from  the  non-local  contributions  that  they  tend  to  fall 
under  different  experimental  peaks  and  thus  final  state  correlation 
effects  are  visible  in  the  Si  I^W  Auger  lineshape.  Although  the 
energy  alignment  (both  absolute  and  relative)  is  rather  good,  it 
does  appear  for  both  the  °2s°2p  anc^  °2p°2p  kands  t^iat  the  separa¬ 
tion  between  the  local  and  non-local  contributions  should  be  about 
2  eV  larger.  This  could  easily  be  accounted  for  by  somewhat  larger 
values  of  the  two  center  repulsion  integrals  u(si,  ,  0-  ) . 

jp  zp 
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3.  Summary  and  Conclusion 

In  both  the  0  KW  and  the  Si  Auger  lineshapes  there  is 

some  evidence  that  the  two  center  repulsion  integrals  u  should 

mn 

be  larger.  The  "classical"  approximation  is  the  Mataga-Nishimoto 
(MN) ,  eq.  (13),  but  others  have  recently  been  suggested.  In  a 
recent  study  by  Dewar  and  Thiel  (35),  the  Klopman  approximation  (K ) 


was  compared  with  the  MN  approximation  and  analytically  evaluated 
two  center  integrals.  The  K  integrals  were  found  to  be  in  general 
agreement  with  the  analytical  integrals,  the  MN  integrals  generally 
too  small.  in  the  event  that  the  two  center  overlaps  are  known, 
Huang  (15)  has  proposed  the  approximation  (h) 


Using  the  overlap  integrals  evaluated  for  Si02  by  Gilbert  et.  al.  (37)  , 
we  obtain  the  values  given  in  Table  4  along  with  the  MN  and  K 
integrals. 

All  three  approximations  behave  properly  in  the  limits  Rmn  -* 

0  and  R  -  o &  .  it  is  in  the  intermediate  regions  that  the 
mn 

approximations  differ.  On  the  other  hand,  it  has  been  argued 

(38,39)  that  the  appropriate  long-range  repulsion  should  behave 
2 

as  e  /(kR  )  where  "k"  is  the  dielectric  constant  of  the  solid 
mn 

and  accounts  for  the  background  electron  density  in  solids.  In 


Si02  k  4,  which  would  significantly  decrease  the  Coulomb  repul¬ 
sion.  However,  in  the  short  range  region  of  interest  here  (near¬ 
est  neighbor) ,  the  macroscopic  dielectric  constant  is  probably 
not  appropriate.  The  results  in  Table  4  indicate  the  MN  integrals 
are  smaller  than  the  others  in  the  short  range  region.  Our 
interpretation  of  the  experimental  Auger  lineshapes  suggest  the 
M  or  K  approximations  to  be  the  more  appropriate. 

In  summary,  we  have  in  this  work  elucidated  the  role  of  final 
state  hole-hole  correlation  effects  in  Auger  lineshapes.  In  the 
0  KLL  case,  the  Auger  lineshape  reveals  a  strong  localization  onto 
a  single  Si20  cluster;  indeed,  it  almost  appears  atomic-like. 

In  the  Si  L23W  case,  the  Auger  lineshape  reveals  significant 
contributions  from  both  the  local  (both  holes  on  the  same  Si20 
sub-cluster)  and  non-local  (the  two  holes  on  neighboring  Si20 
sub-clusters)  configurations. 

Perhaps  most  significant  in  this  work  is  the  use  of  the  Auger 
lineshapes  to  obtain  near  quantitative  information  about  the  hole- 
hole  repulsion  magnitudes  and  bandwidths.  In  general,  the  absolute 
energies  of  the  Auger  transitions  indicate  the  magnitude  of  the 
one-center  integrals,  the  relative  energies  and  intensities  of 
the  local  and  non-local  contributions  some  information  about  the 
two  center  repulsion  integrals  and  bandwidths.  This  suggests  that 
the  quantitative  interpretation  of  Auger  lineshapes  might  be  a 
useful  tool  for  obtaining  important  electronic  structure  informa¬ 
tion  in  other  insulators  and  semiconductors  and  in  chemisorbed 
systems. 
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TABLE  I.  Definition 

of  Atomic  Orbitals  and  Hamiltonian 

Matrix  for 

Two  Model  Clusters 

a.  Bonding  and  Antibonding  Bands: 

B1-A1"B2-A2-33 

Atomic  Orbitals:  <pa  , 

A1 

<PA 

A2 

'"HI  =  (hB1+hB 

)/V2, 

2 

•>H2  ~  2)^2 

Hamiltonian  Matrix: 

I^A1> 

l<PHl)  1  ^2) 

I^A2> 

<^A1  | 

ttA 

V  0 

0 

<*Hll 

V 

°H  r 

0 

<^H2  1 

0 

r  °H 

V 

<<PA2  I 

0 

0  V 

aA 

b.  Non-bonding  Bands  on  A: 

w 

A2-B2-A3_B3~A4 

Atomic  Orbitals:  <pa  ,  <p 

A1  i 

<PA 

A4 

Hamiltonian  Matrix: 

l%>  N>a,> 

1*0 

LfJ 


TABLE  2. 


Summary  of  Matrix  Elements  (tight  binding  parameters) 
obtained  in  this  work  and  comparison  with  other 
reported  results. 


Matrix 

Element 

This  work  (eV) 

Tight  Binding  Calc.  (eV) 

Chadi  et  al14  Nucho  et  al^5 

a(Sih) 

-3  ( -4) a 

-4.6b,c  -7.9b 

a(02g) 

24 

19.  lb 

04 

CM 

O 

b 

6  (5,8)a 

6. 3b  6.0b 

v(°2s'°2s) 

0 

.15 

v(°2p'°2p) 

1 

.45 

v(Sih,  Sih) 

2d 

1.6e 

v<Sih'  °2s) 

5 

6.2® 

v(Sih'°2p  } 

6d 

5.4®  5,  2f 

o 

ro 

01 

15® 

u(°2p) 

15* 

U<SV 

ll,*'h 

aValues  in  parentheses  were  used  for  the  Si405  cluster  where 
a(o2p(b))  and  a(o2p(nb))  were  allowed  to  be  different,  i.e., 

5  and  8,  respectively. 

b4.5  eV  added  to  literature  values  to  change  zero  point  energy 
from  the  top  of  the  valence  band  to  the  Fermi  level  (l). 

ca(Sih)  =  (a(Si3g)  +  3a(Si3p)y4  for  sp3  hybrid  (40).  Values  of 
a(Si3g)  and  a(Si3p)  from  ref.  14. 


( continued) 


(Table  2,  continued) 


dln  the  Si403  cluster  calculations,  a  linear  combination  of  the 
sp3  hybrids  h±  were  utilized,  i.e.,  H  -  ^  (h^hj) .  Thus 
V(SiH,SiH)  =  1/2  V(Sih,Sih)  =  1  and  V(SiH,02p)  =  Y^V (Sih,02p)  « 

eV(Sih,Sib)  =  l/4(o.(Si3s)  -  o  (Si3p)  )  (40). 


V<Sih'°2s/p>  =  1/2[V(Si3s,02  )  +  V(Si  ,0  )j 

O  c  0  c0 

where  by  summetry  V(Si3p  >02syp  )  =  0. 


fTwo  different  sets  of  parameters  were  reported,  the  values 
depending  upon  the  region  of  the  DOS  fitted  (25). 


9. 


gu  was  treated  throughout  as  a  parameter  in  the  calculations  to 
gradually  turn  on  the  effects  of  correlation.  However,  our  best 
estimates  are  given  here  as  determined  from  Manns  (18)  integral 
tables  and  estimates  of  the  correlation  or  static  relaxation 
contributions  (l). 


hu(SiH)  =  l/2[u(Sihl)  +  u(Sihl,Sih2)3  =  7. 
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I  - 


TABLE  4 


The  Auger  Transition  Rates  After  Cl 
from  the  Tetrahedral 


A 

Molecular  Term 


°2i°2s 


0,  -0, 

2p  2p 

terms 


°2s°2p 

terms 


\  (2) 


XT,  (2) 


1/3A,  (4) 


1/3T-  (6) 


Transition  Rate 


[  DCa2, 

1A1)c2(t2)Mcpp(1S)  )2 

[  DUptj.  1T2)cs(«1)Cp<i,1)MC9p(1P) 

+D(t2-  lT2>cp2(t2)('«MCpp(lt,)  '2 

'  c»<oi)yai)Mc»p(3p)  i2 
[  c„(t.)Mpnn(Jp)  r 

P  1  cpp 

[  cp<t2)(-‘)Mcpp(ll,>  )2 
[  1/3A1)c3(a1)cs(a')Mesa(1/3S) 

+0<t2t2-  1/3Ai)cp(t2)Cp(t')Mcpp(1/3s>  ]2 

[  D(.lV  1/3T2)ea(a1)cp(tpKci|p(1'3P) 

+D<V2-  1/3Vc.(Vcp(c2)Mc^<1/3p> 

+D(t2t2,  I2>cp<t2)cp<t2)(-6)Mcpp< 

1  cp(t2>'p<t2>(-4)%P<1/3l’)  ]2 

'  Cp(t2)<:p(t2>%p<1/3p>  I" 


(l/3W  l2 


Number  in  parentheses  indicates  number  of  states  having  this  symmetry. 
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TABLE  5. 


Comparison  of  Approximations  for  the  Two-Center  Coulomb 
Repulsion  Integrals  in  SiOa1 


Item 

Si-0 

0-0 

Si-Si 

R  (A) 

mn  v 

1.5 

2.5 

3.0 

Overlap  <<Pm/<Pn>2 

0.26 

0.063 

0.33 

1  7 

Mataga-Nishimoto  (eV) 

5.5 

4.1 

3.3 

Klopman36  (eV) 

7.7 

5.3 

4.4 

Huang15  umn  (eV) 

10.4 

6.3 

6.8 

''"The  assumed  values  of  u 
respectively. 


(SV 


and  u(o  )  were  11  and  15  eV, 
2P 


2From  Gilbert  et.al.  (37). 
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Figure  Captions 


Fig.  1.  (a)  A  comparison  of  the  experimental  (solid  line)  o  KW 

Auger  lineshape  with  the  theoretical  lineshape  determined 
previously  (l).  (b)  The  difference  between  the  experimental 

and  calculated  lineshapes.  The  peal®  at  475  and  495  eV  have 
been  identified  to  be  shake  satellites.  The  peak  at  511  eV 
is  the  subject  of  this  work.  (c)  A  comparison  of  the  experi¬ 
mental  (solid  line)  Si  L23W  Auger  lineshape  with  the 
theoretical  lineshape  determined  previously  (l).  (d)  The 

difference  between  the  experimental  and  calculated  lineshapes. 
The  peaks  at  52  and  70  eV  have  been  attributed  to  shake 
satellites,  but  in  this  work  are  shown  also  to  arise  partially 
from  correlation  effects.  Possible  sources  of  the  peak  at 
84  eV  are  summarized  elsewhere  (1,27,41,42). 

Fig.  2.  Auger  energies  and  relative  intensities  of  the  eigenvalues 
arising  from  an  LCAO-MO-CI  calculation  on  the  Si403  cluster 
for  various  values  of  uQ  and  ugi#  Contributions  from  the 
°2p(b)-°2p(b) '  O^nbJ-O^b)  and  02p(nb)-02p(nb)  bands  are 
indicated  by  the  light  solid,  open,  and  heavy  solid  lines, 
respectively.  Contributions  involving  the  02s  band  are 
ignored  here.  Numbers  indicate  the  intensity  in  arbitrary 
units. 
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calculation  on  an  0^  cluster.  Columns  two  through  four 
contain  the  relative  °2p(hb)  -C>2p(nb)  Auger  energies  and 
intensities  of  the  eigenvalues  arising  from  a  Cl  calculation 


on  the  01/:  cluster  for  various  values  of  one-  and  two-center 
16 

hole-hole  repulsion  integrals  u^  and  u^-  The  total  intensity 
has  been  normalized  to  100.  The  local  DOS  energies  are  equal 
to  16  eV  minus  the  energies  indicated;  the  Auger  energies  are 
equal  to  517  eV  minus  the  energies  indicated.  The  peak  numbers 
refer  to  Table  3.  The  02p(nb)  bandwidth  as  exhibited  in 
column  one  is  ~6  eV. 


Fig.  4.  Comparison  of  results  for  the  02^(nb)  -02p(nb)  Auger  line- 
shaoe.  On  the  lower  left  is  plotted  the  percent  non-local 
intensity  of  the  total  Auger  intensity  vs.  u/r.  On  the 
upper  right  is  plotted  the  difference  in  energy  between  the 
centroids  of  the  local  and  non-local  Auger  contributions  vs. 
u/r.  u  is  the  one-center  hole-hole  repulsion  energy  parameter; 
r  is  the  02p(nb)  bandwidth  as  determined  from  the  spread  of  the 
one-electron  DOS  of  the  corresponding  n  cluster.  Results  are 
shown  for  clusters  containing  three  atoms  (dotted  line) ,  seven 
atoms  (dashed  line),  16  atoms  with  u^s  =  0  (dot-dashed  line),  and 
16  atoms  with  u. .  /  0  (solid  line). 

i] 
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Fig.  5.  Comparison  of  results  for  the  Si  Auger  lineshape. 

In  column  one  are  the  Auger  energies  and  relative  intensities 
resulting  from  the  fold  of  the  one-electron  DOS  local  to  the 
central  Si  atom  as  determined  from  a  LCAO-MO  calculation  on 
an  Si504  cluster.  The  position  of  the  final  state  holes  is 
indicated  at  the  far  left  in  terms  of  the  M.O.'s  as  named  in 
tetrahedral  symmetry.  The  symbols  (s,  p)  in  parentheses 
indicate  the  dominant  central  Si  atom  contribution.  In 
column  two,  hole-hole  repulsion  is  included  without  Cl  mixing 
for  =  10  and  uq  =  15  eV.  The  third  column  gives  results 

with  Cl  mixing.  The  terms  (A,  E,  or  T)  are  given  as 

determined  by  group  theory  in  T^  symmetry;  the  dominant  con¬ 
figuration  mixing  given  in  parentheses  after  the  term 
symbol.  Numbers  given  indicate  intensities  in  arbitrary  units. 

Acknowledgement;  The  author  wishes  to  thank  the  Surface  Chemistry 
Branch  of  the  Naval  Research  Laboratory  for  its  assistance  in  the 
preparation  of  this  manuscript  and  for  valuable  discussions  with 
several  members  of  the  branch. 
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